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Abstract 

Double stranded DNA chain is known to have nontrivial elasticity. We study the effect of this 
elasticity on the denaturation profile of DNA oligomer by constraining one base pair at one end 
of the oligomer to remain in unstretched (or intact) state. The effect of this constraint on the 
denaturation profile of the oligomer has been calculated using the Peyrard-Bishop Hamiltonian. 
The denaturation profile is found to be very different from the free (i.e. without the constraint) 
oligomer. We have also examined how this constraint affects the denaturation profile of the oligomer 
having a segment of defect sites located at different parts of the chain. 
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I. INTRODUCTION 



DNA is one of the most complex and important biomolecule as it is central to all living 
beings. It contains all the informations needed for birth, development, living and probably 
sets the average life. Structurally it is a giant double stranded linear molecule with length 
ranging from 2/x m for simple viruses to 3.5 x 10 7 /zm for more complex organism [1]. How 
such a molecule came into being during the evolution of life and how it acquired the ability 
to store and transmit the genetic informations are still a mystery. The recent progress in 
genome mapping and availability of experimental techniques to study the physical properties 
of a single molecule [2, 3] has, however, made the field very active from both biological and 
physical point of views. 

A DNA molecule is not just a static object but a dynamical system having rather complex 
nature of internal motions [4]. The structural elements such as individual atoms, group of 
atoms (bases, sugar rings, phosphates) fragments of double chain including several base pairs, 
are in constant movement and this movement plays an important role in the functioning of 
the molecule. The solvent in which the molecule is immersed acts as a thermal bath and 
provides energy for different motion. In addition, collision with the molecules of the solution 
which surrounds DNA, local interactions with proteins, drugs or with some other ligands also 
lead to internal motion. These motions are distinguished by activation energies, amplitudes 
and characteristic times. The motions which are of our interest here are opening of base 
pairs, formation of bubbles along the chain and unwinding of helix (denaturation). The 
energy involved for these motions is of the order of 5-20 kcal/mole. These motions are 
activated by increasing temperature, increasing pH of the solvent action of denaturation 
agents etc. The time scale of these motions are of the order of microseconds and is therefore 
generally unobservable in atomistic simulations as these simulations are restricted to time 
scale of nanosecond because of computational cost. 

The nature of thermal denaturation leading to separation of two strands has been investi- 
gated for several decades [5, 6]. Experimentally a sample containing molecules of a specific 
length and sequence is prepared. Then the fraction of bound base pairs as a function of 
temperature, referred to as the melting curve, is measured through light absorption, typ- 
ically at about 260 nm. For heterogeneous DNA, where the sequence contains both AT 
and GC pairs, the melting curve exhibits a multistep behaviour consisting of plateaus with 
different sizes separated by sharp jumps. These jumps have been attributed to the unwind- 
ing of domains characterized by different frequencies of AT and GC pairs. The sharpness 
of the jump in long DNA molecules suggests that the transition from bound to unbound 
is first-order. The understanding of this remarkable one-dimensional cooperative phenom- 
ena in terms of standard statistical mechanics, i.e., a Hamiltonian model with temperature 
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independent parameters is a subject of current interest [7]. 
II. MODEL HAMILTON! AN AND ITS PROPERTIES 

Since the internal motion that is basically responsible for denaturation is the stretching 
of bases from their equilibrium position along the direction of the hydrogen bonds that 
connect the two bases, a DNA molecule can be considered as a quasi one dimensional lattice 
composed of N base pair units. The forces which stabilize the structure are the hydrogen 
bonds between complementary bases on opposite strand and stacking interactions between 
nearest neighbour bases on opposite strands. Each base pair is in one of the two states: 
either open(non hydrogen bonded) or intact (hydrogen bonded). A Hamiltonian model that 
has been found appropriate to include these interactions and describe the displacement of 
bases leading to denaturation is the Peyrard-Bishop model (PB model) [8]. The PB model 
is written as 

N 

i=i 

2 

H(yi, y i+ i) = + ™{Vi, Vi+i) + V(Vi) (2-1) 

where m is the reduced mass of the base pair, denotes the stretching of the hydrogen 
bonds connecting the two bases of the i th pair and 

( d Vi\ 

The on-site potential V{%/i) describes the interaction of the two bases of the i th pair. The 
Morse potential 

V(y t ) = A(e-° w - l) 2 (2.2) 

which is taken to represent the on-site interaction represents not only the hydrogen bonds 
connecting two bases belonging to opposite strands, but also the repulsive interactions of 
the phosphates and the surrounding solvent effects. The flat part at large values of the 
displacement of this potential emulates the tendency of the pair "melt" at high temperatures 
as thermal phonons drive the molecules outside the well and towards the flat portion of the 
potential. 

The stacking interactions are contributed by dipole-dipole interactions, 7r-electron systems, 
London dispersion forces and in water solution, the hydrophobic interactions. These forces 
result in a complex interaction pattern between overlapping base pairs, with minimum energy 
distance close to 3.4 A in the normal DNA double helix. The following anharmonic potential 
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model mimic these features of the stacking energy: 

W(y u y i+1 ) = ^/c[l + pexp(-«(t/i + y i+1 ))](t/i-^ +1 ) 2 (2.3) 

where k is the coupling constant and the second term in the bracket represents the anhar- 
monic term. When due to stretching the hydrogen bonds connecting the bases break, the 
electronic distribution on the bases is modified causing the stacking interactions to decrease. 
This is taken into account by the exponential term in Eq.(2.3). It may be noted that the 
effective coupling constant decreases from k(l + p) to k when either one of the two inter- 
acting base pairs is stretched. This decrease in coupling provides a large entropy in the 
denaturation. The parameter a in Eq.(2.3) defines the "anharmonic range". 



Spring representing stacking energy 



particle of mass "m" 



V(y) 




"external" potential acting 
on each particle 



FIG. 1: 

Schematic of quasi-one dimensional model of DNA 

The model described above can be viewed as a model of one-dimensional monoatomic lattice 
(see Fig. I) with each atom having mass m and nearest neighbour interaction given by 
Eq.(2.3). Furthermore, each atom is subjected to "external" potential given by Eq.(2.1) 
which effect is to confine the chain in the potential well. The melting takes place because of 
competition between the thermal energy which leads to displacement and "external" field 
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and the nearest neighbour interactions which lead to the confinement. The model, therefore, 
represents a one-dimensional system that differs from the usual one-dimensional systems that 
do not show phase transitions. 

The model Hamiltonian of Eq.(2.1) has extensively been used to study the melting profile 
of a very long (N — > oo) homogeneous DNA chain using both statistical mechanical calcu- 
lations and constrained temperature molecular dynamics [9, 10]. Analytical investigation 
of nonlinear dynamics of the model suggests that intrinsic energy localization can initiate 
the denaturation [11]. The model for a long homogeneous chain exhibits a peculiar type 
of first-order transition with finite melting entropy, a discontinuity in the fraction of bound 
pairs and divergent correlation lengths. However, as the value of the stacking parameter a 
increases and the range of the "entropy barrier" becomes shorter than or comparable to the 
range of the Morse potential the transition changes to second order. The crossover is seen 
at a/a = 0.5 [7]. Though the PB model seems capable of explaining the multistep melting 
in a sequence-specific disorder [12], how this disorder affect the nature of transition has yet 
to be understood. In other work the PB model has been used to understand the melting 
profile of short chains [13] and effect of defects on this profile [14] . 

III. DENATURATION PROFILE 

In a given system of DNA chains the average fraction 9 of bonded base pairs can be written 
as 9 = 9 ext 9 int . 9 ext represents the average fraction of the strands forming duplexes (double 
strands), while 9i nt is the average fraction of unbroken bonds in the duplexes. The equilib- 
rium dissociation of the duplexes C 2 to single strand C\ may be represented by the relation 
C 2 # 2Ci. The dissociation equilibrium can be neglected in the case of long chains as 9 ext 
is practically 1 while 9 int and therefore 9 goes to zero. This is because in the case of long 
DNA chains when 9 goes practically from 1 to zero near the denaturation transition, while 
most bonds are disrupted and the DNA has denatured, few bonds still remaining prevent 
the two strands getting apart from each other. It is only at T >> T m (T m being the melting 
temperature at which half of the bonds are broken) there will be real separation. Therefore 
at the transition the double strand is always a single molecule and in calculation based on 
PB model one has to calculate only 9 int (= 9). On the contrary, in the case of short chains 
the processes of single bond disruption and strand dissociation tend to happen in the same 
temperature range, therefore, the computation of 9 ext in addition to 9i nt is essential. 
Unfortunately, at present, we do not have any reliable method for calculating 9 ext . The 
method which has been used is based on the partition function of rigid molecules and ad- 
justable parameters [5, 13] to be determined from experimental data. To avoid this shortcom- 
ing of the theory we in the present article discuss the denaturation profile of oligonucleotides 
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of given sequence with a base pair at one end of the chain held in such a way that it remains 
at its equilibrium seperation(i.e. there is no stretching) at all temperatures. This can be 
done by creating a deep potential well for this base pair or attaching the one end of both 
strand to a substrate. This will be referred to as chain with constraint in order to distinguish 
it form the "free chain" . One of the advantage of having a constraint of this type is obvious; 
the problem of divergence of the partition function for the PB model for short chains no 
more exists. 

The DNA molecule is known to have non trivial elastic properties: When the two strands of 
DNA molecule are pulled apart by applying a force at one end of the chain, a novel phase 
transition is found (in case of infinitely long chain) to take place at which the two strands 
are pulled completely apart [15]. The phase digram plotted in the plane of temperature and 
force reveals the elastic properties of the DNA chain. Our study reported in this paper differs 
from the situation just described in two ways; (i) a short chain of 21 given base sequence is 
considered and (ii) instead of pulling the chain apart the end base pair is constrained to be 
in unstreched or intact position. 

The oligonucleotide which we consider has the following sequence given by: 

5 ' ACGCT AT ACTC ACGTT AAC AGy 

yTGCGATATGAGTGCAATTGTC 5 ' (3.1) 

The denaturation profile of this oligonucleotide has been studied by Campa and Giansanti 
[13] and by us [14]. We take the same parameters as in the previous study. Thus, D AT = 0.05 
eV, Dec = 0.075 eV, a AT = 4.2A" 1 , a GC = 6.9A~\ k = 0.025 eVA~ 2 , p = 2 and a = 
0.35A- 1 . When one end of chain is held fixed at y — distance, the fraction of intact bonds 
is found using the relation 

1 N 

* = ^5>fa)-yi)> (3.2) 

1=1 

where ~&{y) is Heaviside step function, N the total number of base pairs and the canonical 
average (•) is evaluated by considering all configuration of the chain with one end fixed. The 
i th bond is considered bound if the value of yi is smaller than a chosen threshold y . For 
y we have taken a value of 2 A. Since the model Hamiltonian of Eq.(2.1) couples only the 
nearest neighbours the calculation of canonical average reduces to multiplication of finite 
matrices. The discretization of the coordinate variables and introduction of a proper cut-off 
on the matrices and the number of base pairs in the chain the number of matrices to be 
multiplied. 

In Fig. 2 we compare the value of 6 as a function of temperature for the free and the con- 
strained chain. The influence of the constraint on the denaturation profile is enormous. 
While in case of free chain the breaking of bonds takes place in a narrow temperature range, 
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FIG. 2: 

Variation of 6 and d9 jdT as a function of temperature for free (solid line) and constrained chain 

(dashed line) 



in the constrained chain it is over a wide temperature range. The value of temperature 
at which dd/dT is maximum shifts about 45°C and the peak is much wider and smaller 
compared to that of free chain. 

Next we calculated the denaturation profile for the following two nucleotides having a seg- 
ment of chain with defect sites: 

(a) 5> AC GCTATACTC AC GTTAAC AG 3 > 

yTCGCT AT ACTCTGC AATTGTC^' (3.3) 
(h) 5> AC GCTATACTC AC GTTAAC AG 3 , 

yTGCGATTACTCACGTTTTGTC 5 ' (3.4) 

While both oligonucleotide have ten defect sites their locations differ. In (a) the defects are 
on the left end from site 2 to 11 while in (b) it is in middle from site 6 to 16. The position 
of the base pair is counted from left. The purpose is to see how these defects affect the 
denaturation profile and the formation of loop and stem as often seen in a single strand 
DNA or RNA. In Fig. 3 we plot the variation of 6 and dd/dT as a function of temperature. 
While the denaturation takes place at low temperature the shift is more when the defects 
are in the middle. 

The other quantities of interest are the mean value of displacement of n th base pairs defined 
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FIG. 3: 

Fig. (a): Variation of 6 as a function of temperature for the constrained chain without defect 
(solid line) with 10 defect at one end of the chain (dashed line) and 10 defects in the middle 
(dot-dashed line). Fig.(b): Plot of dO/dT for the same. 



as 

N \ N 



(Vn) y 



\ I \tl d yi)yn e ^p[-P^2 H (yi,yi+i)} ( 3 - 5 ) 
- 1 J \i=i J i=i 



and its fluctuations 

(l^nl 2 ) = \\ (n^l (Vn- (l/n)) 2 exp[-/? £ Hfa, Vi+1 )] (3.6) 



,1=1 / 1=1 



Here Z = J (jjf =1 dy^j exp[— (3 Yl!i=\ H(Vh Vi+i)] i s the partition function of the chain. 
Because of the constraint that the first base pair of the chain is in unstreched condition 
the value of (y n ) as well as (\5y n \ 2 ) depend on the site n. In Fig. 4 we plot the value of 
(y n ) as a function of n at several temperatures for all three constrained chain as described 
above. We find in all cases the opening of the chain starts from the open end of the chain. 
In Fig. 5 we plot (\5y n \ 2 ) as a function of temperature for several base pairs. The quantity 
(\5y n \ 2 ) measures the transverse correlation length for the base pair n. This correlation 
length remains almost zero for all n when the oligonucleotide is in the native state but at 
the denaturation its value increases. At a given temperature the value of (l^nl 2 ) depends 
on n and it increases with n. For a long chain we expect it to diverge for large value of n. 
In conclusion we wish to emphasize that the effect of constraining a base pair at one end of 
a given oligonucleotide has very significant effect on the denaturation profile. 
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FIG. 4: 



Plot of (y n ) vs. n (site position) at four different temperatures. Circles denotes the constrained 
chain without defect, square denotes the chain with 10 defect at one end while diamond denotes 

the chain with 10 defect in the middle. 
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FIG. 5: 

Plot showing the variation of transverse correlation length with temperature for different value of 

n. 
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